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Abstract 

In this paper we discuss possible signatures for the production of scalar t%t* (stopo- 
nium) bound states ai at hadron colliders, where t% is the lighter scalar top eigenstate. 
We first study the decay of <rr ; explicit expressions are given for all potentially impor- 
tant decay modes. If t\ has unsuppressed two-body decays, they will always overwhelm 
the annihilation decays of crj . Among the latter, we find that usually either the gg 
or hh final state dominates, depending on the size of the off-diagonal entry of the 
stop mass matrix; h is the lighter neutral scalar Higgs boson of the minimal super- 
symmetric model. If m<j ( - happens to be close to the mass of one of the neutral scalar 
Higgs bosons, QQ final states dominate (Q = b or t). W + W~ and ZZ final states 
are subdominant. We argue that cr^ —* 77 decays offer the best signal for stoponium 
production at hadron colliders. The tevatron should be able to close the light stop 
window left open by LEP searches, but its mass reach is limited to m a ~ < 90 GeV. 
In contrast, at the LHC one should ultimately be able to probe the region m CTf < 700 
GeV, if the hh partial width is not too large. We also comment on the feasibility of 
searching for production at hadron colliders in the ZZ, Z7 and t + t + t~t~ final 
states, and briefly mention production at 77 colliders. 
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1) Introduction 



Although the Standard Model (SM) of Particle Physics Q has so far passed all experimental 
tests, it has long been known p| to be technically unnatural: Nothing protects the mass of 
the Higgs boson, and hence the scale of electroweak symmetry breaking, from large (quadrat- 
ically divergent) radiative corrections which "naturally" push it up to the Planck scale or 
the scale M x of Grand Unified Theories (GUTs). The probably most elegant solution || of 
this problem is the introduction of N = 1 supersymmetry (SUSY) [|J. In supersymmetric 
theories corrections to the mass of the Higgs boson from loops involving SM particles are au- 
tomatically cancelled by loops involving their superpartners. The cancellation is not perfect 
since supersymmetry has to be broken; naturalness arguments then suggest that the scale 
of SUSY breaking should not (much) exceed 1 TeV. 

So far searches for the direct production of superparticles (sparticles) in collider exper- 
iments have not been successful. However, the lower limits on their masses that can be 
inferred from these searches are only around 120 GeV for most strongly interacting spar- 
ticles 0, and under 50 GeV for all other sparticles [[J. This leaves a wide region to be 
explored by present and future experiments, and much work in that direction has already 
been done @, [jjj. 

Meanwhile various indirect (loop) effects due to supersymmetric particles have been in- 
vestigated. A by now well-known result |§ is that the introduction of supersymmetry allows 
for a beautiful unification of all three gauge couplings of the SM at scale M x — 10 16 GeV. 
In contrast, nonsupersymmetric theories can be unified only at the cost of the somewhat 
ad hoc introduction of new degrees of freedom and/or intermediate scales between Mx and 
My/. Unfortunately these analyses only tell us || [J that the scale of SUSY breaking should 
not exceed 10 TeV or so, and hence offer no immediate clues where to look for more direct 
evidence for the existence of supersymmetric particles. 

Such a clue might come from the third main motivation for the introduction of SUSY. In 
addition to technical naturalness and simple Grand Unification, supersymmetric theories also 
offer the possibility to understand (as opposed to parametrize) electroweak gauge symmetry 
breaking, in terms of (logarithmic) radiative corrections to scalar masses. Even if all these 
masses are identical at some very high energy scale where SUSY breaking becomes effective, 
as e.g. in minimal Supergravity theories 0, radiative corrections will drive the square of the 
mass of one Higgs boson doublet to negative values at low energies, leaving all other squared 
scalar masses positive fl(|. The driving force in this radiative symmetry breaking is the large 
Yukawa coupling of the top quark: Radiative corrections due to Yukawa couplings reduce 
scalar masses, while gauge interactions increase them. This mechanism not only establishes a 
causal link between the breaking of supersymmetry and the breakdown of electroweak gauge 
symmetry, it also points towards a fundamental role for a large (top) Yukawa coupling, and 
might thus eventually help to understand why the top quark is so much heavier than all 
other SM fermions. 

In more practical terms, these considerations indicate that the superpartners of the top 
quark, the stop squarks, might be considerably lighter than the other squarks |IT| ; recall 



that radiative corrections due to Yukawa interactions reduce scalar masses. In addition a 
large Yukawa coupling implies large mass mixing between the superpartners of left- and 
right-handed top quarks, which further reduces the mass of the lighter stop eigenstate £i. 
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As a result, even in minimal supergravity models m t ~ can be almost arbitrarily light even if 



all other squarks have masses of several hundred GeV |L2|, [13 . 

What is the experimental situation? Stop squarks are color triplets, and thus have 
substantial pair production cross sections at hadron colliders. However, present experimental 
bounds |J on squark masses assume 10 or 12 degenerate squark eigenstates, and therefore do 
not apply to ti, at present, searches for events with large missing transverse energy cannot 
exclude the existence of a light stop if the mass of the lightest super symmetric particle (LSP) 
exceeds 12 GeV or so Stop squarks also have electromagnetic and weak interactions. 



However, for a certain mixing angle 9t the Zt\t\ coupling vanishes flql , so that no 

strict lower bound on can be derived from the study of Z decays. In this scenario t\ pair 
production at e + e~ colliders can only proceed via photon exchange. The lower bound on 
mjj will then again depend on the LSP mass, since for too small a t\ - LSP mass difference 
t\ pair events look more and more like hadronic two-photon events.^ 

One might worry that a light stop squark produces unacceptably large loop effects. In- 
deed, a light stop can lead to large corrections to the electroweak p parameter and related 



quantities plfl . However, the main effect would be a reduction of the value of the top quark 



mass mt fitted from electroweak data; note that at present the central value of the SM fit of 
m t |22|] is substantially above the direct search limit on m t [23]. Moreover, this search limit 



is only valid in the SM. If t\ is light enough, the decay t — > t\ + LSP can significantly dilute 
the di-lepton signal for top, especially if m t < 90 GeV In any case, a small need not 



imply large contributions to dp ||21|| . Similarly, loops involving stop squarks and charginos 
can contribute significantly to the matrix element for b — > 37 decays However, this 



contribution can be cancelled by other non-SM loop contributions involving charged Higgs 
bosons. As a result, the recent bounds [f25| on the branching ratio for b — > 57 do not exclude 
|18|| a very light t\. 

We thus conclude that a t\ of 40 or 50 GeV could quite easily have escaped detection so 
far; if we are willing to finetune the Zt\t\ coupling and the ti - LSP mass difference even a 15 
or 20 GeV t\ is not excluded. This also indicates that it is difficult to obtain stringent bounds 
on m% from searches for open t\ production. In particular, most signals will disappear in the 
limit where t\ becomes (almost) degenerate with the LSP; unlike for, say, first generation 
squarks and gluinos this is possible even in the restrictive class of minimal supergravity 
models @, 0.| 



On the other hand, if the t\ - LSP mass difference is small, t\ will be rather long-lived 
27]. The reason is that in this case two-body tree-level decays like ti — > t + LSP or 
ti — > 6 + chargino are kinematically forbidden. In this situation t\ decays preferably into 
a charm quark and a neutralino via a loop diagram, whose decay width is suppressed by a 



*The TOPAZ collaboration at the TRISTAN collider recently observed ]l§ a slight excess of D* mesons 



when trying to measure two-photon production of charm. This has caused speculations 17, 13| |18[ that there 
might be a stop squark with mass between 15 and 20 GeV; tit* production would then give a similar signal as 
two-photon production of charm if the LSP lies just a few GeV below the stop. However, the QCD prediction 
for a (77 — > cc) is quite uncertain fljfl ; the measured cross section is only about 1.5 standard deviations 
above the upper range of QCD predictions. Moreover, preliminary data from the VENUS collaboration 
pil indicate that the ti-LSP mass difference would have to be as small as 2.5 GeV for t\ to have escaped 
detection. Interpretation of the excess in terms of a rather contrived SUSY model therefore seems premature. 

^At least in models with exact R— parity i\ cannot actually be the LSP, since searches for exotic isotopes 
imply fl26|] that the LSP has to be electrically and color neutral. 
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factor ~ 10~ 7 compared to tree-level decays [27 1. Being long-lived, stops can form bound 
states ("stoponia"), which eventually decay via t\t\ annihilation into final states that only 
contain SM particles. 

In this paper we study the decay of scalar (S-wave) t\t\ bound states as well as 
possible signals for production at hadron (super) colliders. Scalar stoponium has been 
studied previously in refs. |2"E|, |2T| |3(J. However, in refs.[^SJ and |JIJ mixing between the 
superpartners of left- and right-handed top quarks was ignored, and diagrams involving 
Higgs bosons in the intermediate or final state were treated only in an approximate fashion 
or not at all; both effects can be very important. They have been included in ref.|29|, but 
there only a very light 01 was treated, so that many decay channels were kinematically 
forbidden. We computed the decay widths for all tree-level two-body decays of cr^ for 
general stop mixing and the whole range of masses of interest in the foreseeable future. As 
already pointed out in our recent Letter j31| the two-photon decay of cr^ probably offers the 
best signal at hadron colliders. Here we present a more detailed discussion of the region of 
parameter space where this signal is viable. We also point out that an interesting range of 
err masses should already be accessible at the Fermilab Tevatron. 

The rest of this article is organized as follows. In sec. 2 we discuss decays. Final states 
consisting of two gluons or two scalar Higgs bosons usually are dominant among decays that 
proceed by tit* annihilation. Contrary to claims in refs. [[2~8| , |30| final states containing two 
heavy gauge bosons contribute much less to the total width of ; this is a direct consequence 
of electroweak gauge invariance. In sec. 3 we discuss signals for production at present 
and future hadron colliders. We focus on the clean two-photon signal, whose background can 
be calculated relatively reliably. At future supercolliders a signal consisting of 4 r leptons 
may also be viable, but here realistic background estimates are much more difficult. Sec. 4 
contains a summary of our main results and some conclusions. Explicit expressions for all 
er^ two-body decays are collected in the Appendix. 



2) Stoponium decays 

In this section we discuss the decays of scalar t\t\ bound states . We work within the 
framework of the minimal supersymmetric standard model (MSSM) which is a super- 
symmetrization of the SM with minimal Higgs sector. In particular, every quark gets two 
superpartners described by complex scalar fields, which are associated with the left- and 
right-handed components of the quark Dirac spinor. If the quark is massive, its two super- 
partners can mix with each other. For the case of stop squarks the resulting mass matrix is 
given by (in the basis (tL,tn)): 

2 / m| + m\ + 0.35m| cos2/3 — m t (A t + /i cot/3) \ 

* = \ L -m t {A t + fj, cot/3) m\ R + m\ + 0.16m| cos2/3 ) ' ( 1 ' 



where we have used the conventions of ref.|[L2|, which are quite similar to those of ref.[]3"2| 



The quantities m^, mi R describe the soft SUSY breaking contributions to the diagonal 
squark masses. As already mentioned in the Introduction, loop corrections involving the top 
Yukawa coupling tend to reduce these mass parameters compared to their values at higher 
energies. In models where all squark masses are equal at some very high (GUT or Planck) 
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scale one therefore expects mi L to be smaller than the corresponding quantities for squarks 
of the first two generations. The same argument also implies < mi L , since Yukawa 
interactions affect the running of mi more strongly [fTOf . 

The diagonal entries of the stop mass matrix (JI|) also depend on the mass m t of the top 
quark as well as the angle (3, defined via tan/3 = (H%)/(Hf), where Hi, H 2 are the two 
Higgs doublet fields necessary in any realistic supersymmetric model ||. Two additional 
parameters enter the off-diagonal entries of ([]]): The supersymmetric Higgs(ino) mass fi, 
and the parameter A t describing the strength of the trilinear scalar H^t^t^ interaction, 
which breaks supersymmetry. A t and \x are free parameters of the model, but we generally 
expect them to be of roughly the same magnitude as mi h ; tan/3 can be taken to be positive, 
but both At and \x can have either sign. 

The mass matrix (|l|) can easily be diagonalized. We define 




cos6* t sm8 t 
— sin^t cos^t 




and obtain 



11,2 



m LL + m RR =F V ( m LL ~ m RR) + ^™LR 



(2) 



(3a) 



m\ - m LL 



taxi9 t = — 1 — 2 > (3b) 



m LR 

where iti 2 lllrrr are the (1,2), (1,2) and (2,2) elements of the mass matrix (|1|), respectively. 
Note that cos 2 t < (>) 1/2 if mf > (<) mf . 

We will see below that cr^ decays can involve charginos and neutralinos as well as Higgs 
bosons; we therefore briefly describe the corresponding sectors of the MSSM. The charginos 
Wi, W 2 are Q mixtures of the superpartners of the W ± bosons and of the charged Higgs 
fields contained in Hi 2 . Similarly, the neutralinos Z i: i = 1,...,4 are mixtures of the 
superpartners of the (unmixed) B and W 3 gauge bosons of the SM as well as of the neutral 
Higgs bosons. Charginos are Dirac fields while neutralinos can be described by Majorana 
spinors. The mass matrices for charginos and neutralinos depend 0] on the parameters 
(jl and tan/3 introduced above, as well as on the SUSY breaking gaugino masses Mi and 
M 2 . We will for simplicity assume the usual unification relation Mi = 5/3 tan 2 8\yM 2 — 
M2/2; similarly, the gluino mass is given by M3 = a s sin 2 9wM 2 /a. The description of 
the neutralino, chargino and gluino sectors thus only necessitates the introduction of one 
additional parameter, which we take to be the mass M 2 (at the weak scale). Finally, at 
the tree level the Higgs sector of the MSSM is described by two parameters ||32|| , which 
we take to be tan/3 and the mass mp of the pseudoscalar Higgs boson. We have included 
radiative corrections to the Higgs masses and mixings from top-stop loops |33[], employing 



the effective potential formalism [34]; these corrections are determined completely in terms 
of the parameters entering the stop mass matrix (|I]). 

We are now in a position to discuss <ji decays. There are two very different kinds of 
decay modes: Single stop decays, where one stop squark decays independently of the second 
squark; and annihilation decays, which proceed by tit\ annihilation. The partial widths 
due to single stop decays are simply twice the corresponding stop decay widths. Since we 
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are interested in relatively light stop squarks, we will assume that t% — > t + g decays are 
kinematically forbidden. We then have to consider the following two-body t\ decays: 

fa -> bW h i = 1,2; (4a) 
fi->t%j = l,...,4; (4b) 
fi ->c2T i ,j = 1,...,4. (4c) 

The processes Q4"a|,b) occur at tree level, with full gauge or top Yukawa strength. If they 
are kinematically unsuppressed, the corresponding t\ decay widths are very roughly of order 
10~ 3 to 10 -2 77if ; exact expressions are listed in the Appendix, eqs.( |A.16| , [A.17f ). 



If the channels (|4a|,b) are kinematically closed, fl4"c|) is usually the dominant decay mode 
27fl of t\. Note that it couples a third generation squark to a second generation quark, via 
a neutral (super) current. In models where the squark mass matrix and the quark Yukawa 
coupling matrix commute at some energy scale, e.g. if all squarks are mass degenerate 
at some energy scale, the process flic] ) therefore has to proceed via a weak loop diagram 
involving small elements of the quark flavor mixing matrix. Even though the amplitude is 



enhanced by a factor \og(Mx/Mw) the authors of ref.p7| therefore estimate the squared 
matrix element for the process (|4c[) to be suppressed by a factor ~ 1CT 7 compared to the 
tree-level process Q4"b"D, leading to a decay width ~ lO^m^. We will see below that this 



is completely negligible compared to the widths for annihilation decays, to which we turn 
next. 

As already mentioned, these proceed via the annihilation of the t\ and t* that make 
up erf ; this kind of decay by far dominates the total widths of the familiar lowest lying 
quarkonium states (r) c , J/tp, T). Here we are only interested in two-body decays of erf , 
which dominate all other annihilation decays. We treat nonrelativistic bound state, 

where the squarks are in an S— wave. The partial width for er^ — > A + B is then given by 



3§: 



3/5 -l^(o)l 2 — - 

Here, 



Tfa^A + B) = ^J^R^^-^- E \M^»m^AB)\l =0 . (5) 



\ 



m\ + m 2 B \ Am\m 2 B 
m 2 ) mi 



(6) 



is the usual phase space factor, and 1/(1 + 5ab) is a statistics factor. Finally, R(0) is the 
wave function at the origin. For realistic QCD potentials the wave function generally has to 
be computed numerically. Ref.[^7f provides parametrizations for R(0) as well as the binding 



energy of the first ten S— wave states of a nonrelativistic (s) quarkonium system, using a 
potential that describes the known cc and bb systems well. We will use their parametrizations 
throughout. 

Eq.(|5J) reduces the problem of computing ct£ annihilation decay widths to the calculation 
of the Feynman amplitude A4 XaXb for the annihilation of t\ + t\ into A + B with helicities 
and Xb, respectively. Here, the initial state is assumed to be a color singlet, and summation 
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over color degrees of freedom of the final state is understood.]] Since we are only interested in 
S— wave initial states we need the Feynman amplitude only in the limit of vanishing relative 
velocity v of the stop squarks; this simplifies the calculation considerably. 
We computed the matrix elements for the following processes: 





-> QQ\ 


(7a) 


n~ — 


_> w + w~- 


I ' "J 


a h ~ 


-> ZZ; 


(7c) 


°tl " 


+ Z T ; 


(7d) 


a h - 


-> 77; 


(7e) 


a h - 


-> hh; 


(7f) 


a h - 


-> 66; 


(7g) 


°h ~ 


-> tt; 


(7h) 


°h ~ 


-> ZiZj, i, j — 1, . . . , 4, 


(7i) 



where h in eq.([7D stands for the light neutral scalar Higgs boson. In general t\tl annihila- 
tion can proceed via the four classes (topologies) of Feyman diagrams shown in fig. 1. The 
t— channel diagram la contributes to all modes of eq.(^), but with different particles being 
exchanged: t\ for the gg, Z'j and 77 final states; t\ and t 2 for the ZZ and hh final states; 
for the W + W~ final state; a chargino for the bb final state; a neutralino or gluino for the tt 
final state; and a top quark for the Z^Zj final state. The u— channel diagram of fig. lb only 
contributes if the final state particles do not carry any conserved charge (reactions |7a],c- 
f,i); the exchanged particles are then the same as for the t— channel diagram. The 4-point 
diagram of fig.lc only contributes if the final state particles are bosons (processes |7a|-f). 
Finally, in the limit v — > only the neutral scalar Higgs bosons h, H can be exchanged 
in the s— channel diagram of fig. Id; this diagram therefore only contributes to reactions 
(fTE|,c,f-i). Explicit expressions for the matrix elements for the processes of eq.(0) are listed 
in the Appendix. 

A first example of cr^ branching ratios is shown in fig. 2. For clarity not all final states 
listed in eqs.(ffl), (0) are represented in this figure. We have fixed m t = —fi = 150 GeV, 
m ih = m iR = 200 GeV, tan/5 = 2, m P = 500 GeV and M 2 = 100 GeV, and have varied A t 
between -67 and -312 GeV; since A t and \x have the same sign, m a ~ decreases monotonically 
with increasing \A t \, see eq. (|3~a|) f\ 

We see that for this choice of parameters the by far dominant decay mode of is into 
two gluons, as long as the single stop decays (f£a|,b) are kinematically forbidden. In fig. 2 the 
mass of the lighter chargino is around 110 GeV. For m CT - > 2 (rn^ + m^j ~ 230 GeV the 

* Following ref.|36| the color wave function of the initial state in eq.(||) has been normalized such that 
the color factor is 1 if A and B are singlets under SU(3); this explains the factor of 3 in eq. @. See the 
Appendix for more details. 

Un this and the following figures we have set m a - i — 2mi i , i.e. have neglected the small reduction of m ai 
due to the binding energy. This is consistent since the treatment of refs.]35|, |3(| also neglects the binding 
energy when computing cr^ decay widths. 
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decay ([la]) (not shown) opens up and quickly dominates over all annihilation decay modes 
(0). Notice, however, that even for m„. = 400 GeV the width of the single stop decay mode 
([4~aD is "only" a few hundred times larger than that for the gg final state; since the partial 
width for the decay mode ([4c]) is ~ 10~ 7 smaller than that for ( fja|) we conclude that this 
loop decay is entirely negligible as far as decays are concerned. 

Given the large widths for the single stop decay modes Q4"a|,b), one might worry whether 
our treatment is adequate for these decays. Inherent to our approach is the assumption 
that formation and decay of can be treated separately. This is only true if the at 
formation time is significantly shorter than its lifetime. A good measure for the (inverse 
of) the formation time is the binding energy Eund of ■ For a purely Coulombic potential 
the time required to complete one Bohr orbit is proportional to l/Eund, and we expect this 
relation to survive qualitatively also for the more realistic QCD potential of ref.[|37|]. In fig.3 
we therefore compare the binding energy of the lowest (IS 1 ) cr^ state, as parametrized in 
ref. f37j, with the total decay width, for two different choices of parameters. In both cases 
we took m t = 150 GeV, tan/3 = 2, m h = 400 GeV, m~ tR = 300 GeV and m P = 500 GeV; 
77V- was varied by changing A t . The solid curve has been obtained with /x = —300 GeV and 
fixed M.2 = 100 GeV. In this case the light chargino is mostly an SU(2) gaugino (wino), and 
the lightest neutralino is mostly a U(l) gaugino (bino). Recall that mi > mi R implies that 
ti is dominantly an SU(2) singlet (cos 2 9 t < 1/2), so that in this scenario the tiWib coupling 
oc cos9 t is suppressed; cos 2 9 t decreases with decreasing iz, — £r mixing and increasing in 
this case, which explains the flattening of the solid curve at large m a ~. We see that in this 
case the total cr^ decay width is still a factor 3 to 10 below the binding energy, even well 
beyond the threshold for decays (|4a|,b). Our results of fig. 2, where W\ is also mostly a wino, 
should therefore be at least approximately correct. 

In contrast, the dashed curve has been obtained for fixed /i = —80 GeV, while M 2 has 
been increased along with mj . The lightest chargino and neutralino states are now both 
higgsino-like, so they couple to II and Ir with equal (top Yukawa) strength. The total cr^ 
decay width therefore increases rapidly with m a ~\ moreover, the opening of the t% — > t + Z\ 
channel is more pronounced than in the previous case. As a result, r(cr t - i ) does indeed become 
comparable to the binding energy in this scenario, which means that our approach will not 
work for m a - > 300 GeV or so. Methods that have recently been developed to describe the 



ti threshold f38[ will have to be adopted instead. 

Finally, in both cases we observe a very prominent peak at m a . = 505 GeV, where the 
s— channel heavy Higgs boson exchange diagrams become resonant. Since the total decay 
width of the heavy Higgs boson exceeds the binding energy, a proper description of this 
case would have to combine the methods of refs.]38|] with the results of refs. |[40|| where the 
mixing between a nonrelativistic bound state with a (narrow) Higgs resonance is discussed. 

The results of fig. 2 show that the branching ratios for those annihilation decays that 
might yield a detectable signal for er^ production at hadron colliders (see sec. 3) become 
very small if the single stop decays ([|a|,b) are allowed; fig.3 showed that cr^ may not have 
time to form at all if there are light higgsino-like states. In fig. 4 we have therefore chosen 
our parameters such that these single stop decays are kinematically forbidden for m a - < 600 
GeV; this has been obtained by choosing [i = —300 GeV, with the other parameters having 
the same values as for the dashed curve in fig.3. Comparing fig. 4 with fig. 2, we notice two 
obvious differences. One is the structure around m a - = 505 GeV, which is due to H exchange 
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becoming resonant as already discussed in connection with fig. 3. Of course, in the immediate 
vicinity of the resonance our results are not reliable, but the ti final state remains dominant 
in regions of parameter space where r(crf 1 ) is well below the binding energy. 

The other prominent feature of fig. 4 is the very rapid increase of the branching ratios 
for the hh and, to a lesser extent, W + W~ and ZZ final states. In case of the hh final state 
this can be explained from the observation that in the relevant limit m 2 p 3> m 2 z the ht\t\ 
coupling |41] contains a term cx m t (A t + \i cot/3) /mw oc m 2 LR /mw, where m 2 LR is again the 
off-diagonal entry of the stop mass matrix ([I]). Moreover, in fig. 4 mi L and mi R are rather 
large, so that the hh threshold occurs at a point where m 2 ^ <C m| ; eq.(|3a]) shows that this 

also implies m 2 LR ~ min(m 2 ,m~ R ) 3> m^. Close to threshold the t\ exchange contribution 
therefore scales like 

rnin(m~ , m~ ) 

M(htl -> hh)\ h exchange oc L R m\ > m 2 , m\ (8) 



see eq. (|A.9|) ; the square of this amplitude clearly decreases very quickly as we move away 
from the hh threshold (oc w^ 4 )- This rapid rise of r(cr t - i — > hh) at threshold has first been 
observed in ref . p9| ,f] 

The behavior of the matrix elements for the W + W~ and ZZ final states is somewhat 
more complicated. In the region m| ^> m\y we can use the equivalence theorem |^2| 



to understand the production of longitudinal gauge bosons, which in this region usually 
dominates the production of transverse gauge bosons. This theorem states that Feynman 
amplitudes involving external longitudinal gauge bosons are the same (up to corrections of 
order l/m^y) as those where the gauge bosons are replaced by the corresponding would-be 
Goldstone bosons (GBs). These Goldstone modes can also be produced via the Feynman 
diagrams of fig.l. However, it is important to note that the squarks which are exchanged in 
the t— and u— channels are heavy here. In case of the neutral GBs only t% contributes since 
there are no diagonal couplings between pseudoscalar Higgs bosons and stop squarks J32 



The corresponding diagrams for charged GBs contain a bi squarkp]; note that m~ bL is linked 
to mi L by gauge invariance: 

m\ L = m 2 L - m 2 w cos2/3. (9) 

The nonvanishing trilinear scalar couplings of the GBs have similar strength as the ht\t\ 
coupling; the corresponding t— and u— channel diagrams are therefore at best of order 

min(?77/~ , Tn~ 1 

MitJl -> GBs)\ q exchange OC % , m\ > 111 2 , m\ , (10) 

i.e. are suppressed by a factor of order (m^/m^ compared to those for the hh final state. 
In addition, there are s— channel h exchange contributions 

_ min(?77/~ , Tfi~ ) 

MiUl -> hh, GBs)\ h exchange OC - 2 ^ t R , m\ > m\, TT& (11) 

4m| — m^ L ' R 1 



■fin that paper the ti exchange contribution to hh production has not been included; this contribution is 
small where the hh channel is important. 

^Diagrams with 1>r exchange are |32| proportional to ml and can thus be neglected. 



This contribution exists for both hh and W + W and ZZ final states, but is suppressed by 



a factor of order (mw l m t L r) compared to the t\ exchange contribution (||) to hh produc- 
tion. Far above threshold the width for longitudinal gauge boson production is therefore 

suppressed by a factor of order (rn^/m^ compared to the width for the hh final state. 

The equivalence theorem is not applicable close to the W + W~ and ZZ thresholds. Near 
the thresholds the t— and u— channel diagrams, which only contribute for longitudinal gauge 
bosons as shown in eq. (|A.3|) , are suppressed by powers of the phase space factor (3 of eq.(f|). 
The 4-point and s— channel h exchange diagrams therefore dominate here; the curves of fig.2 
show that they often interfere destructively. Note that for m| <C m~^ , eq.flllf) also applies 
approximately for the s— channel contribution to W + W~ and ZZ production; it is this term 
which leads to the rapid increase of the corresponding partial widths near threshold.^] 

Why did refs. |28| and [30 find so large branching ratios for the W + W~ final state? 
The crucial omission is that the relation (|9[) between mi and mj £ has not been taken 
into account in these papers. We emphasize again that this relation follows directly from 
SU(2) gauge invariance (and its spontaneous breakdown); it is independent of the details 
of supersymmetry breaking. A violation of eq.(^) therefore implies explicit (hard) gauge 
symmetry breaking, which renders the theory non-unitary and/or non-renormalizable. In 
addition, t^-tR mixing has been neglected in these papers. Looking at the exact expression 
flOp for M (titj -> W+W-) it is clear that the cancellation between the t— channel and 
4-point diagrams which ensure unitarity in the simple limit cos 2 9 t = 1, = m^ L = mi 
can be spoiled if one choses m~ ^> . However, this implies either cos 2 6t — > (if m 2 LR 
is kept fixed), or a large s— channel h— exchange contribution which again "conspires" to 
restore unitarity; recall that gauge invariance relates the ht\t\ coupling to the - t R mass 
splitting and hence to if mj £ ~ m t R - I n short, the suppression of the partial widths for 
a t 1 W + W~, ZZ is a textbook example for the unitarity restoring cancellations that are 
so characteristic for gauge theories. 

Figs. 2 and 4 show that the branching ratios for all other modes listed in eqs.fl7|) are 
quite small. The width for the 77 final state is simply 8a 2 /(9a 2 ) ■ r(crf l — ► gg), and the 
partial width for the Z7 final state is of similar magnitude or even smaller. (Recall that 
the Zt\t\ coupling vanishes JOJ for cos 2 6 t = 4/3 sin 2 9w) For the parameter choices of 
figs. 2 and 4 the light neutralinos are gaugino-like; the partial widths for the ZiZj final 
states are therefore comparable to those for the 77 and Z7 final states. Had we chosen 
the light neutralinos to be higgsino-like, their partial widths would have been larger by a 
factor oc {m t /mw) A - For parameter choices leading to mixed-state neutralinos (where both 
gaugino and higgsino components are sizable) the s— channel h exchange contribution to 
ZiZj production can become important, leading to partial widths comparable to those of 
the W + W~ and ZZ final states. However, the existence of light higgsino-like or mixed-state 
neutralinos in the MSSM also implies a small mass for the light chargino, so that the single 
stop decay mode ( 4a ) is allowed, totally swamping all ox annihilation decay modes as we 
have seen above. 

Finally, the partial width for the bb final state is very small unless m a . ~ or tan/3 3> 1. 
At first glance the gaugino exchange diagram seems to contribute with full SU(2) gauge 



^The coupling of the heavy scalar Higgs boson H to W and Z bosons is suppressed Q for trip ^> m 2 



W ■ 
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strength. However, chirality implies that r(cr t - 1 — > //) oc m 2 for any SM fermion /. As a 
result the bb final state can be important only if the b— quark Yukawa coupling is enhanced 
(tan/? ^> 1), or in the immediate vicinity of the h resonance. Fig. 4 shows that even the 
partial width for the ti final state is quite small away from the H pole. This is partly due 
to destructive interference between s— and t— channel diagrams for m a , < mjj, and partly 
because color factors suppress all t— channel contributions, see eq.( |A.14| ). 

We thus conclude that, if the tree-level single stop decays fl4a|,b) are kinematically for- 
bidden and m CTf is not close to either m h or m#, the total a% decay width is dominated 
either by the gg or hh partial width, with the W + W~ and ZZ partial widths playing a 
secondary role. Our discussion of fig.4 already showed that the ratio of the gg to hh partial 
widths crucially depends on the size of the LR element of the stop mass matrix (|]). This is 
further illustrated in fig. 5, where we show the gg, hh and W + W~ partial widths as a func- 
tion of m a ~ for three different choices of parameters. We have fixed /i = 500 GeV, mp = 1 
TeV and m t = 150 GeV; instead of varying A t we have fixed the dimensionless quantity 
A = 2 A t /{mi + J7i£ ) as well as the ratio m^/m^, and varied mf . The dependence of 
the hh and W + W~ partial widths on m a ~ is therefore quite different than in figs. 2 and 4|] 
In particular, eq.(|5a) implies that now m 2 LR , and hence the strength of the ht\t\ coupling, 
increases with increasing m a . even if <C m f LR - F° r large A the hh partial width still 

decreases with increasing m av due the t\ propagator suppression, but the decline is much 
less rapid than in fig.4. Notice that we have chosen fi > here, so that m 2 LR increases 
monotonically with increasing A. However, due to destructive interference between different 
diagrams both the hh and the W + W~ partial width initially decrease with increasing m 2 LR , 
shooting up quickly once A > 1. We have already seen above that the partial width for the 
W + W~ final state always stays well below those for the hh and gg final states. Here we see 
that the off-diagonal entries of the stop mass matrix (|l|) need not be all that large for the 
hh mode to dominate oi decays. Finally, the short dashed curve has been obtained with 
mi L = rrit R , as compared to mi R = 0.7 mi h for the other curves; we see that this has only 
little effect on T{pi x — > hh). Since we kept A fixed, the size of m 2 LR for given is about 
the same for the two choices of mi R / mi L ; this again indicates that the size of m 2 LR is indeed 
the quantity that decides whether or not the hh partial width is sizable. We will come back 
to this point later. 

Having completed our discussion of decays, we are now ready to study possible signals 
for its production at hadron colliders. 



3) Signals for stoponium production 

In this section we discuss how one might search for cr^ states. We focus on hadron (su- 
per) colliders, since they offer the largest cross sections; furthermore, the machines we dis- 
cuss either already exist (tevatron) or are in a relatively advanced stage of planning (LHC), 
while plans for future linear e + e~ or 77 colliders are still at a very preliminary stage. The 
production of cr^ at hadron colliders proceeds via gluon fusion. This process is related by 
crossing to the cr^ — > gg decay, whose partial width we computed in the previous section. 

"Of course, the gg partial width is fixed uniquely by m CT - and the strength of the QCD coupling constant, 
and his hence the same for all three cases. 
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The total cross section for o~ tl production is then (to leading order in the QCD coupling 
constant) simply given by 

a(pp -> a h + X) = I' -tG(x, Q 2 )G (-, Q 2 ) ^T(a h -> gg), (12) 
Jt x \x J am%, 

where r = m 2 a Js with s being the squared pp invariant energy; since in leading order only 
gluon fusion contributes, the cross section is the same for pp and pp colliders. Unless stated 



otherwise we have used the parametrization of ref.|43| for the gluon distribution functions G, 
and have chosen the scale Q = m$ . In order to set the stage for the subsequent discussion we 
remark here that eq . fll2|) predicts a total oi x production cross section at the LHC (\/s = 14 
TeV) of 80 (22, 8.5) pb for m CT _=150 (200, 250) GeV. 

Unfortunately most ai decays will not lead to a signal that is detectable at hadron col- 
liders. To begin with, the QCD di-jet cross section integrated over any reasonable invariant 
mass window will be many orders of magnitude larger than the total production cross 
section, making it impossible to detect — ► gg decays. QCD backgrounds also swamp 
<7f — ► bb, ti decays. Single stop t% — > b + W\ decays can give hard, isolated leptons in the 



final state if W\ decays leptonically. However, the open t\t\ pair production cross section [13 
is some 4 orders of magnitude larger than the ai cross section. The presence of 2 invisible 
LSPs and at least one neutrino in the — ► W^W^bb — > / ± X signal makes it impossible 
to reconstruct m a ~ even in principle. We conclude that single stop decays will never give a 
signal for cr^ production at hadron colliders even if conditions are favorable for the detection 
of open stop production |L4 



In ref. |HJ the use of the W + W~ decay mode (for rather heavy a^, to be produced at 
future supercolliders) was advocated. However, we have seen in the previous section that in 
this paper the Br{ai 1 — > Vy + Vy~) was overestimated by a large factor. Moreover, this final 
state can only compete with QCD backgrounds if both W bosons decay leptonically. The 
event will then contain two neutrinos, which make it impossible to reconstruct the invariant 
mass of the W + W~ system. Even if it were possible to somehow discriminate against the 
enormous ti background at hadron supercolliders, the continuum cross section for W + W~ 



production is still at least an order of magnitude larger [44] than the ax — > W + W signal. 
The authors of ref.[p9| proposed to use the cr^ — >• hh — > t + t + t~t~ decay as a signal. 



This might have been feasible for light and light h at the tevatron collider. Unfortunately 



the lower bound on nih from Higgs searches at LEP j|5] excludes this possibility for light 
enough to be produced at the tevatron. The t + t + t~t~ SM background is much larger at 
supercollider energies, making it considerably more difficult to observe a cr^ signal in this 
channel. For example, at the LHC (y/s = 14 TeV) the ZZ — > 4r background amounts 
44] to approximately 10 fb. Using [46| Br(h — > t + t~) = 10% we find that the err — > 4r 



signal could be as large as 1 pb if m r < 150 GeV and Br(az i — > hh) ~ 1. While this is 
considerably larger than the most narrowly defined physics background, it is smaller than 
the cross section for pp — > bbbb —>■ At production^. Backgrounds from events containing 
jets with low charged particle multiplicity, which might fake r signals, could be even more 
dangerous. Note that the presence of (at least four) neutrinos in the final state makes it 

*This is true for the high luminosity option of the LHC, where most Ab events originate from independent 
pp collisions. The 46 — * 4r background for low luminosities isj47] around 0.2 pb. 
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once again impossible to reconstruct either rrih or m a -. Isolation of the t + t + t~t~ signal at 
a hadron supercollider therefore seems quite difficult to us; certainly detailed Monte Carlo 
simulations would have to be performed before a good case for this signal can be made. 

This leaves us with ai decay modes with rather small branching ratios. The Z^Zj mode 
again suffers from the problem that m a . cannot be reconstructed since the final state con- 
tains two LSPs. The ZZ mode offers a clean signal if both Z bosons decay leptonically. 
Unfortunately the branching ratio for the ZZ final state is often quite small, as shown in 
figs. 2 and 4. The conditions for this signal are most favorable for large LR element of the 
stop mass matrix (p]) and in the mass range rrih > > mz- In this case Br{ai 1 — > ZZ) 
can be as large as 10%, giving a maximal a(pp — > cr^ — > l + l + l~l~X) ~ 8/6 for m„. = 200 
GeV at yfs = 14 TeV, corresponding to 800 events in a full LHC year (/ Cdt = 100 fb -1 ); 
this should be readily detectable. However, since in the MSSM rrih cannot be larger than 
140 GeV or so even after the inclusion of radiative corrections p3| , |34|j this window of op- 
portunity is rather narrow. Figs. 2 and 4 show that more generically the ZZ branching 
ratio lies between 0.1 and 1%, making this signal rather marginal; recall that in the SM, 



a(pp -> ZZ -> 1+l+l-rX) ~ 40 fb @ at the LHC. 

The Z'j final state could also give a clean signal if Z — > l + l~ . Unfortunately the combined 
branching ratio for o~ tx — > Z^ — > / + / _ 7 is always below 0.01%. Note that photons have to be 
quite energetic to yield a potential signal at the LHC; this final state can therefore only be 
used for m a - t > 150 GeV or so. The total pp — > — > Zj — > Z + Z~7 signal then amounts to 



at most 10 fb at ^fs = 14 TeV; in comparison, the SM physics background |48[ is about 600 
fb even if one requires the transverse momentum of the photon to exceed 50 GeV. It seems 
therefore very unlikely to us that this signal will be detectable. 



Such considerations led us to propose [311] the 77 final state as the most promising 



signal for cr^ production at hadron colliders. Figs. 2 and 4 show that the corresponding 
branching ratio is typically a few times 10~ 3 , although it can be substantially smaller near 
an s— channel pole or for large rn 2 LR ; this is considerably larger than typical branching ratios 
into leptonically decaying Z bosons. Since the — ► 77 partial width is determined uniquely 
by mi (for given QCD potential), the 77 signal rate depends on model parameters only via 
the total o"^ decay width. The signal is very simple, consisting of two hard photons with 
invariant mass M 77 = m a ~ in a hadronically quiet event. Of course, there is also a sizable 
SM background from qq annihilation and gg fusion. It has been studied in some detail in the 
literature [f4£], |5(| as a background to a possible signal for intermediate mass Higgs boson 



production. Recall that the natural width of the signal peak in our case is just a few MeV, 
see fig. 5; in contrast, the background gives a smooth distribution in M 77 . The question is 
then if, or under what circumstances, the signal peak is observable on top of the background. 

In most SM 77 events the photons will emerge at small angles, due to t— and u— channel 
quark propagator effects; in contrast the signal is isotropic in cos^*, where 9* is the scattering 
angle in the 77 center-of-mass system. We therefore impose the cut 

|cos0*|<O.5. (13) 

The 77 background has been computed to next-to-leading order in QCD [|50|| . However, if 
one vetoes against the presence of hard, central jets in the event and requires the photons 
to be isolated, the NLO prediction for the background rate is actually very similar to the 
leading order estimate. Moreover, no NLO calculation for the signal cross section exists as 



12 



yet. We therefore also treat the background in leading order, but we include the gg — > 77 
contribution which has been found to be very important |49j especially for low M 7 



'77 ' 



As noted above, the natural width of the signal peak is extremely small; however, due 
to detector resolution effects its actual (measured) width will be much larger. Clearly the 
background should be integrated over this larger range of M 77 . On the other hand, the signal 
within a given bin need not be larger than the background in order to be detectable, since the 
expected background level can be determined experimentally by fitting a smooth function 
to the sidebins. The question is then whether the excess in the signal bin is statistically 



significant or not. Following ref.^TJ, where the search for SUSY Higgs bosons was discussed, 
we define the signal to be significant if the 99% c.l. upper limit on the background rate is 
smaller than the 99% c.l. lower limit on signal plus background combined. In the limit of 
large event numbers, where Gaussian statistics can be used, this means: 



N b + N s - 2.32^/ N b + N s >N b + 2.32y/N b , (14) 

which implies 

N s > 2.32 (2^ + 2.32) . (15) 

Here, N b and N s are the expected number of signal and background events after cuts. As 
noted earlier, the background has to be integrated over a detector-dependent bin width 



AM 77 : 



N h = Cdt 



do b 



M 77 =m 



AM 77 . (16) 



In the limit \fW b ^> 1 the minimal detectable signal cross section cr™ n = N s / (J Cdt) therefore 
scales like the inverse square root of the integrated luminosity and also like the inverse square 
root of the energy resolution of the electromagnetic calorimeter, which determines the size of 
AM 77 . In our background estimates we simply took AM 77 to be twice the assumed invariant 
mass resolution. 

In fig. 6 we show the expected signal at the tevatron (y/s = 1.8 TeV). In addition to the 
cut (|T3"| ) we have required that both photons have rapidity \y y \ < 1.1, so that | cos6* 7 | < 0.8 in 
the lab frame; the same cut has been applied by the CDF collaboration in their preliminary 
analysis |)2] of events with two hard photons. The dashed curve has been obtained under the 
assumption that the total width of is determined by the gg and 77 partial widths alone, 
while for the solid curve all cr^ decay modes of eq.fl7[) have been included. The two results 
are indistinguishable except for m a ~ ~ rrih = 87 GeV for the given choice of parameters. The 
signal for such light does therefore not depend on the details of the (s)particle spectrum 
(aside from mj ,) as long as t\ has no tree-level two-body decays and s— channel h— exchange 
contributions to decays are not "accidentally" enhanced. 

In fig. 6 we also show our estimates for the minimal detectable signal cross section (dotted 
lines) for three different values of the integrated luminosity representing the present status 
(18 pb" 1 ), the hoped-for luminosity after run lb (100 pb _1 ), and an estimate of what might 
be achievable after the new Main Injector has been completed (1 fb _1 ). Here we have 
assumed an invariant mass resolution of 2%, i.e. AM 77 = 0.04M 77 . Since the expected 
number of background events per bin is not always large we have used Poissonian statistics 
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to derive these curves; however, eq.fll5|) gives quite similar results in the mass range where 
detection of might be possible. 

We conclude from fig. 6 that the mass reach of the tevatron for err searches is quite 
modest. It is therefore exceedingly unlikely that one of the two CDF events J52[ with very 
large M 77 (350 and 430 GeV, respectively) is due to production. Indeed, most of the 
region that one might be able to probe even with an integrated luminosity of 1 fb _1 is already 
excluded by LEP measurements of the total Z decay width and/or squark searches at LEP, 
unless the Zt±t* coupling happens to be suppressed by — t^ mixing. On the other hand, 
fig. 6 also reveals that existing tevatron data might already help to probe this light stop 
window; in particular, they might confirm or rule out the recently proposed explanation ]l7| 
of the (small) excess of events containing low— D* mesons observed |16| by the TOPAZ 
collaboration in terms of t\t\ production and subsequent ti — > c + Z\ decay with ~ 15 
GeV and m 2l ^ 12.5 GeV. 

In fig. 7 we show results for LHC energy (y/s = 14 TeV). We have again applied the cut 
( |T3"D on cos9*, but have relaxed the requirement for y 7 slightly compared to fig. 6; our cut 
1 2/7 1 — 1-74 still ensures that the photons are at least 20° away from the beam pipes, i.e. 
are well isolated from the beam remnant jets. We have also rather optimistically assumed a 
1% resolution for the measurement of M 77 , i.e. used eq.([T6|) with AM 77 = 0.02M 77 . Since 
now the expected number of background events per bin is quite large we have used Gaussian 
statistics to estimate the minimal detectable signal cr™ n ; the dotted curve shown in fig. 7 
assumes one nominal LHC year of operations, i.e. J Cdt = 100 fb . 

The solid line in fig. 7 shows the total production cross section (divided by 100), 
without any cuts. The dashed curves show the 77 signal cross sections after cuts for the two 
sets of parameters chosen in figs. 2 and 4. We saw in sec. 2 that the branching ratio for the 
77 mode is about 3 • 10~ 3 at small m CTt -, where the gg mode dominates. Comparing the solid 
and dashed lines in fig. 7 we see that our cuts reduce the signal by approximately a factor of 
3.5 at low m U{ ; at high vn a - almost half of all — > 77 events pass. Of course, the cut ( |i~3D 
alone excludes 50% of all signal events; for large m a . the events are automatically central 
so that the cut on the rapidity does not reduce the event number further. The reduction of 
the qq, gg — > 77 backgrounds by our cuts is much larger; in addition we have to require the 
photons to be well isolated from all jet activity, including the beam remnant jets, in order 
to suppress the bremsstrahlung background, which otherwise dominates [fSfl . 

The short dashed curve shows that, as anticipated, the 77 signal quickly becomes un- 
observable once i\ has tree-level two-body decays (m ff . > 230 GeV; see fig. 2). Here the 
situation is quite analogous to the case of ti bound state production where the 77 signal also 
becomes inaccessible [|53| once m t > 120 GeV or so. The rapid decline of the long dashed 
curve in fig.7 at m ai = 180 GeV is due to the opening of the hh mode, see fig. 4. In this case 
the 77 signal becomes marginal just beyond the hh threshold, but should still be observable 
after several years of LHC operations if our assumptions about the detector resolution can 
be realized. Note that the signal actually increases with increasing m a ~ as we leave the hh 
threshold region, in spite of the rapid decrease of the total cross section for production; 
this once again illustrates the steep decline of the hh partial width with decreasing m 2 LR , 
which corresponds to increasing in this case as discussed in sec. 2. Finally, for m a , ~ m# 
the 77 signal again becomes unobservable, due to the large s— channel enhancement of the 
ti partial width. 
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Before we try to further evaluate the stoponium discovery potential of the LHC it might 
be worthwhile to discuss some of the uncertainties inherent to our calculation of signal rates. 
As stated earlier, the cross sections shown in figs. 6 and 7 have been computed using the 
parametrization of the parton distribution functions in the proton given by Owens [43 



we found that other recent parametrizations |54] give very similar results. Our cut on the 
rapidity of the photon ensures that we probe the gluon density G at comparatively large 
values of x where differences between existing parametrizations are not so large. In our 
previous figures we have taken Q = mf for the momentum scale in the gluon distribution 
function; the same choice has been used for the solid and dashed curves in fig. 8. In contrast, 
the dotted curve has been obtained with with Q = m a „. Clearly this change of Q by a 
factor of 2 has little effect on the predicted signal rate. For small (large) values of x the 
gluon density increases (decreases) as Q is increased; most of our signal comes from the 
cross-over region in between, where G depends very weakly on Q. We note here that we 
have not changed the scale in a s in fig. 8, which determines T(a^ — > gg) and hence the total 
a$ production cross section, see eq.fll^p. However, for the given choice of parameters the 
total ui decay width is always dominated by the gg partial width, so that the branching 
ratio for the 77 final state is inversely proportional to r(ay — > gg); the signal is therefore 
almost independent of the choice of the momentum scale to be used in a s in the given case. 

Nevertheless our results do depend on the choice of the QCD scale A, as also shown in 
fig. 8. The reason is that larger values of A imply a bigger QCD coupling constant a s , and 
hence a more tightly bound stoponium system, i.e. larger |i?(0)| 2 ; note that the signal is 
oc |i?(0)| 2 if the total decay width is dominated by annihilation decays. Ref.pTfl provides 
parametrizations of this quantity for four different values of A; our previous results have 
been obtained with A = 0.2 GeV, which is in between the extreme choices of 0.1 and 0.4 
GeV. We see that even for r tot (cr t - i ) ~ r(cf — > gg) the variation of A corresponds to a 30% 
uncertainty of our signal. This uncertainty is even larger if T tot (a^) 3> T(a^ 1 — > gg). If the 
total width is dominated by annihilation decays into hh or ti final states, the uncertainty in 
A leads to an approximately 50% uncertainty of the signal, since now the increase of the gg 
partial width, i.e. of the total cross section for production, is no longer cancelled by a 
corresponding decrease of the branching ratio for the 77 final state when A is increased. The 
A dependence becomes stronger yet if tree-level two-body decays of t% are possible, since 
in this case the signal is oc T(a^ 1 — ► gg)Y{ai 1 — > 77)/T tot (cr t - i ) oc a 2 |i?(0)| 4 ; the signal now 
increases by more than a factor of four when A is increased from 0.1 to 0.4 GeV. A similarly 
strong dependence on A was observed in ref. [53| for the analogous case of the 77 signal 
for toponium production. However, fig. 7 shows that detection of err at the LHC becomes 
much more difficult if r tot (cr^ 1 ) T(ai — > gg), and all but impossible if t x has unsuppressed 
tree-level two-body decays. The situation depicted in fig. 8 is therefore more characteristic 
for situations where the discovery of seems feasible at the LHC. 

Yet another uncertainty comes from the existence of higher (excited) stoponium states. 
So far, we have only considered the direct production of the lowest lying (n = 1) state. 
However, already for the cc system two (J = 1) S— wave bound states are known to exist; 
there are three J = 1 s— wave bb bound states. For a Coulomb potential the number of 
bound states increases proportional to the square root of the mass of the heavy (s) quark. As 
mentioned earlier, in ref. ||37[| the mass (binding energy) and wave function at the origin of 
the first ten heavy (s)quarkonium states are given (the ground state and nine excited states). 
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Not all of these states will be true bound states; some may be resonances that decay rapidly 
into a pair of stop-flavored (spin-1/2) "mesons". We assume, rather conservatively, that the 
mass of the lightest stop "meson" lies just 200 MeV above m^, i.e. we treat the n— th t\t\ 
state oi x (n) as a bound state if 

m a . {n) < 2m h + 0.4 GeV. (17) 

Of course, the wave function at the origin, and hence the production cross section, is 
smaller for excited states than for the ground state (|i?(0)| 2 oc 1/n 2 for a Coulomb potential). 
Nevertheless fig. 9 shows that the enhancement of the signal due to the production of excited 
states can be quite substantial. In this figure we show J2 n l-^n(0)| 2 /|-Ri(0)| 2 , where the sum 
runs over all states that satisfy the condition (|T7|), i.e. are true bound states. This ratio is 
equal to the enhancement of the signal due to the production of excited stoponium states if 
annihilation decays dominate the total decay width, and if the excited states themselves 
contribute to the signal in the same way as the ground state does. 

This second requirement deserves a brief discussion. Since all excited states in the sum are 
S— wave states, expression (|5|) also describes their annihilation decays (with -R(O) replaced 
by R n (0)). The various stoponium states should lie close enough together that their 77 
invariant mass peaks will not be distinguished from each other by the detector, unless we 
have underestimated the resolution to be achieved by future experiments. Annihilation 
decays of the excited states will therefore contribute to the signal in the same was as for the 
ground state. However, these excited states have additional decay channels: They can decay 
into lower lying stoponium states, plus a photon or a mesonic system with vanishing charge 
and strong isospin. Cascade decays of excited states into lower lying S— wave states will also 
contribute to the signal if the lower lying state decays into two photons, since the existence of 
additional very soft photons or mesons from the cascade will hardly be detectable at hadron 
colliders. 

However, an excited state can also decay into a lower lying state with different angular 
momentum, e.g. into a P— wave state. The relative branching ratios for annihilation decays 
of these higher spin states can differ significantly from those of the S— wave states; eq.([|) 
is no longer applicable here. In particular, P— wave states do not contribute to the signal 
at all (unless they in turn decay into another S 1 — wave state), since they cannot decay into 
two photons. In order to estimate how much the higher S— wave states contribute to the 
signal one would therefore have to follow all their decay chains; this necessitates a complete 
understanding of stoponium spectroscopy, which is well beyond the scope of this paper. 
We mention here that only a rather small fraction of T(2s) and T(3s) mesons decays into 
P— wave bb (xb) states which do not decay back into S— wave states (~ 13% for T(2s) and 
~ 22% for T(3s)). If this result carries over to the stoponium system, fig. 9 should give a 
good estimate for the enhancement of the signal due to the production of excited states. 
Finally, there is a contribution to the signal from the direct production of states with J ^ 
and their subsequent decay into S— wave states, but it should be quite small.Q 

^In the most simple treatment of nonrelativistic bound states^, |36|, the production cross section of 
higher spin states is predicted to be small since it is proportional to the square of derivatives of the wave 
function at the origin, divided by additional powers of . Recently Bodwin et al. ifsE) have suggested that 
the production of P— wave quarkonia states might be enhanced by the presence of a sizable component of the 
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In any case, in the absence of a more reliable treatment of the decays of excited squarko- 
nium states we have conservatively decided to only include the direct production of the 
lowest (n = 1) S— wave state in our estimates of signal cross sections. 

We finally address the question of the o~ tl discovery potential of the LHC, using the 77 
decay mode and our conservative estimate of the signal cross section. We have already stated 
repeatedly that will be unobservable at hadron colliders, and may indeed not form at 
all, if the single stop decay modes (f|a|,b) are unsuppressed or cr^ is very close in mass to 
one of the two scalar Higgs bosons of the MSSM; for the subsequent discussion we therefore 
assume that this is not the case. We saw in sec. 2 that the size of the 77 branching ratio 
is then almost uniquely determined by the partial width for the hh final state, the ratio of 
the gg and 77 partial widths being fixed by QCD. Moreover, we saw in fig. 5 that, at least 
for parameters where T(ai — > hh) > r(crj 1 — > gg), the partial width for the hh final state 
is determined by the size of the LR element of the stop mass matrix ([I]). Under the given 
assumptions the detectability of the cr^ — > 77 signal at the LHC therefore basically depends 
on two parameters: The mass m CT ~, which determines the total erj production cross section; 
and m 2 LR = —m t (A t + /j, cot/3), which determines the size of the branching ratio of the 77 
decay mode. 

It can safely be assumed that the mass of the top quark will be known quite precisely 
before LHC experiments are ready to search for o£ production. In fig. 10 we therefore show 
the region in the plane of m ai and A t + /i cot/3 that can be probed after one and five nominal 
LHC years (C = 100 fb _1 per year). Here we have assumed m t = 150 GeV and tan/3 = 2, 
but this choice has little affect on the accessible region.^ 

The general shape of the curves is easy to understand. At small m„. the hh mode is only 
open if A t + fi cot/? is very large; note that radiative corrections reduce if A t + \i cot/3 ^> 



rriT 



33], p4fl . In this case the hh partial width is very large just beyond the threshold, as 
shown in fig. 4, and the 77 signal remains unobservable even after a long running period. On 
the other hand, for larger m cr{ the hh channel is always open. We see that the curves for the 
maximal accessible A t + fi cot/? become quite flat in this region. The reason is that increasing 
m a . decreases the total cross section for oi x production, but at the same time decreases the 
branching ratio for the hh mode if A t + fi cot/3 is kept fixed, see eq.(f|). Moreover, the 
minimal detectable signal cross section decreases with increasing m a ~, although more slowly 
than the total production cross section does, as shown in fig. 7. These effects tend 
to cancel each other, leading to the observed flattening of the curves for m„. > 220 GeV. 
Eventually, however, increasing m cr{ reduces the hh partial width to a value below the gg 
partial width; decreasing it even further then has little effect on the signal, and the curves 
terminate rather abruptly. We finally note that the little bulge in the accessible regions at 
A t + /i cot/? ~ 600 GeV occurs because for moderate values of m 2 LR the hh partial width no 
longer grows monotonically with A t + fi cot/3, as we already saw in fig. 5. 

wave function where the QQ is in a color octet state. However, the same component would also suppress the 
branching ratio for P— wave states to decay into S— wave states, so that altogether the P— wave contribution 
to the 77 signal is still small. 

*It does affect the size of the region in the top left corner where the LEP Higgs bound EH] is violated. 
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4) Summary and Conclusions 



In this paper we have studied the decays of S— wave tit* bound states , as well as possible 
signals for their production at hadron colliders. We first argued in sec. 1 that there are no 
strict bounds on which hold both for all ih — t R mixing angles 9t and all values of the LSP 
mass; even under relatively mild assumptions a ti as light as 40 GeV is still allowed. This 
leaves a wide mass region to be explored. We have seen that production is only detectable 
at hadron colliders if ti has no unsuppressed tree-level two-body decays. Otherwise single 
squark decays of dominate over annihilation decays, and stoponium production gives at 
best a small contribution to the signals for open stop production. 

The dominant annihilation decay modes of oi x are those into two gluons, two light scalar 
Higgs bosons h, or a ti pair. Since the latter two decays involve electroweak rather than 
strong couplings, their partial widths have to be enhanced dynamically in order to be com- 
parable to or larger than the one for the gg final state. In case of the ti mode this can only 
happen if m CT . is very close to the mass of the heavy scalar Higgs boson H , so that s— channel 
H— exchange contributions become (almost) resonant. The hh partial width becomes large if 
the off-diagonal entry m 2 LR of the stop mass matrix is approximately as large as the diagonal 
entries of that matrix; in such a situation mixing greatly reduces the mass of the lighter stop 
eigenstate. Since the hiii\ coupling increases with m 2 LR while decreases, thereby further 
enhancing t\ exchange diagrams, the hh partial width is very sensitive to ttl^ R) as illustrated 
in fig. 5. 

Unfortunately we saw in sec. 3 that none of these three potentially dominant final states 
leads to a readily detectable signal at hadron colliders. The most promising mode appears to 
be the oi — * 77 decay, which gives rise to a peak in the two-photon invariant mass spectrum. 
We analyzed this signal in some detail, comparing it to the 77 continuum background. We 
found that existing Tevatron data might already begin to close the light stop window left 
by LEP (where the Zt\t\ coupling is suppressed by mixing and the ti-LSP mass difference 
is small). On the other hand, even fo r / Cdt = 1 fb _1 the mass reach of the tevatron only 
extends to m„~ = 90 GeV. Under favorable circumstances this mass reach can be extended to 
500 (700) GeV after one (five) year(s) of operation at the LHC with full luminosity (C = 100 
fb _1 per year). Recall, however, that for m a . > 120 GeV the hh decay mode of might be 
open, which might greatly reduce the branching ratio for the 77 final state. More generally 
LHC experiments will therefore only be able to probe a region in the (m ap m 2 LR ) plane, see 
fig.10. 

We should remind the reader here that our calculation has considerable uncertainties, 
even beyond those intrinsic to any leading order QCD prediction for hadronic processes. 
On the one hand, we have ignored backgrounds from jets with very few charged particles, 
which could fake a single photon. This background is clearly detector dependent, but could 
potentially be sizable. On the other hand, our estimate for the signal rate is probably also 
too low, since we have ignored all contributions involving higher (excited) stoponium states. 
We saw in fig. 9 that they might enhance the signal by as much as a factor of two; however, a 
quantitative treatment of their contribution requires a detailed understanding of the entire 
stoponium system. 

Once (j^ production has been observed in the 77 channel its mass will be known precisely. 
If 160 GeV < m ai < 300 GeV one might then be able to find evidence of production also 
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using the 7Z and/or ZZ channel, where Z bosons decay into e + e~ or fi + fi~ pairs; at least 
this task should be easier than searching for production in these channels before m a ~ is 
known. Once the existence of o~z has been established one might also try to look for its hh 
decay via the t + t + t~t~ final state. We saw that the cross section for this final state could 
be as large as 1 pb; the main problem here is to cleanly identify the r leptons. Data taken 
at lower luminosity are probably more useful for this purpose, since the presence of multiple 
overlapping events will make r identification even more difficult. 

Once m a - is known, one can even contemplate studying it in some detail at a 77 collider. 
At least in principle such a device can be constructed |56| by backscattering laser photons 
off the electrons and positrons of an e + e~ collider. The cross section for cr^ production 
could be of the order of (0.5 pb)/(m cr -/100 GeV) 3 . Moreover, by polarizing the incident 
photons one can greatly reduce backgrounds; e.g., 77 — > qq production would be suppressed 
for light quarks if both photons have the same polarization, which might even allow to detect 
a ii ~* 99 decays. The strong dependence of many partial widths on model parameters (see 
figs. 2 and 4) makes their measurement either at a pp or a 77 collider very interesting, and in 
particular offers one of the few possibilities to measure the size of the trilinear soft breaking 
parameter A t . 

Notice that searches for stoponium production are in some sense complementary to 
searches for open stop production. Stoponium states will be very difficult to detect, and 
might not form at all, if t\ decays via two-body modes that are accessible at tree level. On 
the other hand, open stop production at hadron colliders will be difficult to detect either via 
its semi-leptonic decay or via a missing p? signal unless the ti-LSP mass difference is siz- 
able. These two requirements are complementary because within the minimal SUSY model 
there is a strong correlation between the possibility of tree-level two-body decays of ti and 
a large ti-LSP mass difference. This is obvious for the stop — > top + LSP decay, but also 
holds if the t\ — > b + chargino decay is allowed, at least in the case where the LSP (which we 
always assume to be the lightest neutralino) is dominantly a gaugino. A gaugino-like LSP 
is favored dynamically in models with radiative gauge symmetry breaking [M, as well as 
by cosmological considerations; unlike a higgsino-like or mixed-state LSP, it can naturally 
explain the observed Dark Matter in the Universe [57j. 

We thus conclude that there should be a sizable ti-LSP mass difference, which facilitates 
detection of open stop production, if t± decays rapidly; if the ti-LSP mass difference is small, 
the light stop is usually long lived and chances for stoponium production should be good. 
This complementarity is not perfect. On the one hand, the possibility of a large branching 
ratio of az into hh or, worse, ti final states means that we cannot derive a firm "no-loose" 
theorem for stop searches at hadron or e + e~ colliders. On the other hand, if ti is rather light 
it might well be long-lived even if the stop-LSP mass difference is large, since for chargino 
masses below 100 GeV or so the rule of thumb that the chargino is twice as heavy as the LSP 
(for gaugino-like LSP) need not apply. In such a scenario both open stop and stoponium 
production might be observable at hadron colliders, the former via the t\ — > c+LSP loop 
decay, the latter in the two-photon channel. Given the intimate connection between stop 
squarks and the puzzle of electroweak symmetry breaking, in particular in models where this 
breaking occurs radiatively, experimental searches for any signal for scalar top production 
at present and future colliders are well worth the effort. 
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Appendix 

In this appendix we list the (squared) matrix elements for S— wave color singlet tit\ pair 
annihilation into the two-body final states of eqs.^, as well as the widths for the single stop 
decays of eqs.([|a|,b). We do not include annihilation into two charginos, since this can only 
occur if the t\ — > W\ + b decay is allowed, in which case it swamps all annihilation decays. 
Recall also that we assume t\ to be the lightest strongly interacting supersymmetric particle. 

We start with the annihilation decays. In the following expressions we have suppressed 
color indices, i.e. summation over colors has been performed. The resulting color factors 
are included explicitly. As already briefly noted in sec. 2, the color wave function in the 
initial state of the matrix element in eq.([5]) is given by |5o6j after contraction with the 
color indices of the scattering amplitude, a sum over a and b has to be taken. Notice that 
this is not normalized to unity, which explains the appearance of the color factor of 3 in 
eq.fl^); however, this normalization has the practical advantage that annihilation into two 
color-singlet particles has color factor one. 

gg, 77 final state 

For an S— wave initial state, i.e. for v — > 0, these final states receive contributions only from 
four-point interactions (fig. lc); the t— and u— channel t\ exchange diagrams (figs. la,b) 
vanish in this limit (for physical, i.e. transverse, gauge bosons). The squared gg amplitude 
can be written as 

32 

E \M{txt\ -> gg)\ 2 v=0 = —gt- (A.1) 

color, spins 

The squared amplitude for the 77 final state is given by 

E IMMi-77)lLo = 16g 4 e 4 , (A.2) 

spins 

where q = 2/3 is the charge of stop. 
W + W- final state 

This final state receives contributions from diagrams with t— channel exchange of sbottoms 
(only bi contributes if terms cx m& are neglected), s— channel exchange of light (h = H2) 
and heavy (H = Hi) neutral scalar Higgs bosons, and also from the 4 point interaction of 
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two stops and two W bosons. The s— channel exchange of the Z boson does not contribute 
to the S— wave amplitude. Here we list the amplitude for specific helicities A, A of the W 



bosons; A, A can take the values 0, ±1. 



1 2 2 a V- 9H,WWcf 

-fl COS ft - ) K 7T- 



2 



M X \U\ - W+W-) = (SnSiofo + (-I)^aa 

- 2/^7^ao^ 2 cos 2 6 t m , + (A.3) 

The other combinations of helicities do not contribute; this is easily understand from spin 
conservation. In eq.( |A.3|) we have introduced jw = m i 1 / m w and (3w = yl — { m w l™^) 2 '-, 
jw is a kinematical factor which appears in the polarization vector of longitudinal gauge 
bosons (A = 0). 

We have included mixing between SU(2) doublet and singlet stops defined as 

in eq.(Q) of the main text; however, we ignored sbottom mixing. The QHiWW are the Higgs 
W + W~ couplings ||]: 

9HiWW = gm w cos(/3 - a), g H2 ww = gm w sin(/3 - a). (A.4) 



The Higgs tit \ couplings c| are defined in eqs.(A3)-(A5) of ref.|4l[]; we list them here for 
completeness: 

cf = ^^s^( l - cos 2 9 t -l sin 2 9 W cos2ft) + 



ll cosQw 2 3 ' m w 

gm t sin26 t 



u 



2m 



w 



where 



s 



(i) 



:(A t r^+^), (A.5) 



cos(a + /3); = sin(a + /3); 



r (D = sma . r (2) _ cosa . 
sin/3' u sin/3' 

,/n cosa ,/ 9 i sina 
sm/3 sm/3 

Here tan/3 is the ratio of vacuum expectation values introduced in sec. 2, and a is the 
mixing angle of the neutral scalar Higgs bosons @. Note that -> 1 and r^ (2) -> cot/3 
if the pseudoscalar Higgs boson is much heavier than mz] the last term in eq. (|A.5| ) is thus 
proportional to the off-diagonal entry of the stop mass matrix in this limit, as emphasized 
in the text. 



ZZ final state 

The contributing Feynmann diagamms are similar to those for the W + W~ final state, except 
now the t— channel exchanges proceed through tip, and crossed (it— channel) diagrams have 
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to be added since the Z bosons don't carry a charge. We find: 
M X \U\ -> ZZ) V=0 = -^f^ (8 XQ 8- X0 f3 2 z + (-) X 5 X ~ X ) 
1 



cos 2 9w 



(») 

2g T - o 8111 M cos 2 t + - sin 4 W )-Y, T~2 V 

v 4 3 9 / v 4m~ - mi- 



4 
9 



COS z Uw 



(cos 2 9 t - | sin 2 W ) 2 C o S 2 t sin 2 ft. 



2m~ — m| 



*2 



(A.7) 



where ^ z = / m z and (3 Z = \Jl - {m z /m ii y 



Z7 final state 

Since the photon does not have longitudinal polarization states, the t— and u— channel 
exchange of t\ again disappears in the v — > limit. Furthermore, only \ z = A 7 = ±1 states 
are allowed. After summing over the final state polarization, we get 



E \M(tJl - Z 7 )|L, = (^cos 2 ^ - ^sin 2 ^ 

spin COS ^ U 6 



(A. 



hh final state 



Here all four classes of diagrams depicted in fig.l contribute: 



2(4f) 2 



Mm;^«!)U, = j^44: 



m~ + m~ — m 2 ^ 



, r 22 



H =— - — 5 — [2 sin2o; sin(/3 + a) — cos(/3 + a) cos2o;l 

4ra~ - 2 cos 



c- 2 ^ 3 
4 m l- m l 2 2lZ Z w COS2aSinW + a \ 



(A.9) 



Here c~ 2 ~ 2 and c 22 are the ti — 1 2 — h and t\ — t\ — h — h coupling, respectively; they can be 
expressed as 



(2) m z sm(a + (3) . (2 . 2 1 

41 =0 7 sm20 t o sm W - - 

(1 ' 2 costly \3 4 

gmt (a ■ \ oa 
^{A t cos a — fisma) cos 29 t ; 

ny/ sin p 



C ll - o 



2mw sin /3 

cos 2a 1 2 2 2 m? cos 2 a 

— ——(-cos f t --sin 9 w cos29 t ) 5 — r^— 

cos 2 6*vf 2 3 77% sm 2 /3 



(A.lOa) 
(A.lOb) 
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ZiZj final state 

This process proceeds by t— and u— channel exchange of a top quark and s— channel exchange 
of Higgs bosons. In our convention neutralino eigenstates are obtained by diagonalizing the 
neutralino mass matrix by a real orthogonal matrix, thus the mass eigenvalue of a neutralino 
can be either positive or negative. Defining h and h to be the helicities of the two neutralinos 
(h, h = ±1/2), we have: 



M(t x t\ - ZiZj = S h ^iml - {m~ Zi + m^) 2 

2m i (a i aj — bibj) + (rag . + m^.)(ajaj + bibj) 



\{m\ +m|. 



+ 9 



(sin aS'-j — cos aQ'- 
4m^ — m 2 Hi 



mi 



.(2), 



+ 



sin aQ" + cos aS" 
4??!^ — m 2 H2 



(ATI) 



Here a, and 6j are scalar and pseudoscalar stop-top-neutralino couplings; explicit expressions 
are given in eqs.(3), (8) and (9) of ref. JO]. Q^- and S 1 ^- are Higgs-neutralino couplings defined 
in ref. [B2| ; recall that they are real in our notation. 



bb final state 

This process proceeds via the t— channel exchange of charginos (Wi, W 2 ) as well as s— channel 
exchange of Higgs bosons: 



M(htl -> bb) 

2 



„2 



mi 



it 3 m 2 - m| - m 2 ^ 



2mvy cos (5 



(i) 

c~ cos a 
ti 

4m| — m 2 ^ 



+ 



(2) . 

c~ sin a 

£1 

4ra~ — m 2 ^ 



(A.12) 



Here /i and h are again the final state helicities, and q and <ij are scalar and peudoscalar 
stop-bottom-chargino couplings defined as 



hWib 



Ci 



b(ci + da^Wih + h.c; 

g,r a gm b u i2 

--V ll cos6 t + — = - 

l lm w cos p 



di = -~VaCOs9 t 



gm b U t 



1 2 



cos t + 
cos 6 t + 



gm t Vi 



1 2 



2\f2m w sin (5 
gm t V i2 



sin 9 t ] 
sin 9 t . 



(A. 13a) 
(A.13b) 

(A.13c) 



Uij and Vij are 



2\/2mw cos /3 2\/2mw sin f3 

the matrices that diagonalize the chargino mass matrix. We use the 



same conventions as in eqs.(A6) to (A. 8) of Ref. [58]; in particular, rriw l2 can have either 
sign. Notice that c 2 — d 2 is proportional to m?,; hence the whole cross section is suppressed 
by the square of the bottom quark mass. This decay mode therefore turns out to be negligible 
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for the whole parameter space of our interests, unless it is "accidentally" enhanced by an 
s— channel (h or H) pole. The factor of 1/3 for the Wi exchange term is a color factor 
necessary for t— channel color singlet exchange. We have included an overall color factor of 
y/3, which strictly speaking only occurs after the (incoherent) summation over the final state 
color indices. 



tt final state 



This process proceeds t— channel exchange of neutralinos and gluinos (g) and s— channel 
exchange of scalar Higgs bosons: 




-25 hh V3 
-&2) +mt (a? + &? 



ml 



gm t 



2 



2m w sin [3 



2 

- m% 

(1) • 

Cr ' sin a 
ti 



+ 



6|)+m«(af + 6f) 



4m? 

ti 



+ 



(2) 

c} ' cos a 



2 

mf 
ti 



mf 



Hi 



4m; 



m 



#2 



(A.14) 



The couplings a^, fej have already occured in eq. (|A.ll|) above; a§, 6^ are the corresponding 
gluino-stop-top couplings: 



a § + b g = "93 sin26» t . 



(A.15) 



Here again the Z^ exchange term receives a color factor of 1/3, while the g exchange con- 
tribution comes with a factor of 4/9 for color octet exchange. Notice that we have again 
included an overall factor of which properly only appears in the squared amplitude after 
summation over the ti color states. 



Single stop decays: t\ 



W t + b, Zi + t 



If the stop mass is larger than + m^ or m% + m t , single stop decays dominate and the 
pair annihilation modes described above all have a very small branching ratio. The decay 
width into b + W is given by: 



rfo -> bw; 



\M\ 



m~ 



mi 



where 



2c 



in 



+ 2£ 



m 



i - ( m b - 



(A.16) 



(A.17) 



spin 



The couplings c$ and di have been defined in Eq. (|A.13 ). Eqs.( |A.16| , |A.17| ) also describe the 
decay width for t\ — > Zj+t, with the following substitutions: m b — > m t , m-fir. — ► m z^ c « ~^ a j 
and di — > bj\ the tt\Zj couplings aj and bj have already been introduced in eq.( A.ll|) . 
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Figure Captions 



Fig. 1 The four classes of Feynman diagrams that contribute to annihilation decays of stopo- 
nium into two-body final states. 

Fig.2 Branching ratios for annihilation decays of cr^ listed in eq.(|7]). The range of mi x values 
shown results from varying A t between -312 and -67 GeV. The values of the other 
parameters are: = m; = 200 GeV, m t = —fi = 150 GeV, M 2 = 100 GeV, 

mp = 500 GeV and tan/3 = 2. The branching ratios for the bb and ti final states (not 
shown) are always below 10~ 3 . 

Fig. 3 The binding energy of the lowest stoponium state Of (dotted) is compared with the 
total cr^ decay width (solid, dashed), for two different sets of parameters. We have 
chosen m t = 150 GeV, m P = 500 GeV, tan/5 = 2, m h = 400 GeV, and m- tR = 300 
GeV. The solid and dashed curves correspond to scenarios with a gaugino-like and 
higgsino-like LSP, respectively. 

Fig.4 Branching ratios for annihilation decays of cr^ listed in eq.([|). The range of 

values shown results from varying A t between 440 and 1080 GeV. We have increased 
the SU{2) gaugino mass M 2 along with so that the tree-level single stop decays 
of eq.(^J) remain kinematically forbidden (M 2 = 1.5raf ). The values of the other 
parameters are: m t = 150 GeV, mi L = 400 GeV, = —fx = 300 GeV, mp = 500 
GeV, and tan/3 = 2. The branching ratio for the bb mode is again small. 

Fig. 5 Dominant partial widths for ai x annihilation decays. The gg partial width (dotted) 
depends only on m% and the QCD scale parameter A, while the hh (solid, short- 
dashed) and W + W~ (long-dashed) partial widths in general depend on all parameters 
entering the stop mass matrix of eq. ([!]) . We have kept A = 2A t / {yn~ tL + vn- tR ) as well 
as the ratio m^ L /m^ R fixed and varied m^. Most curves are for m^ R = m t ~ L , but the 
short-dashed curve has been obtained with mi = 0.7mr . The values of the other 
parameters are: m t = 150 GeV, tan/3 = 3, // = 500 GeV and m P = 1 TeV. 

Fig. 6 Cross section for pp — > cr^ — > 77 after cuts at the tevatron. The dashed curve assumes 
rtotC^ti) — r(c£ — > gg), while the solid line includes all channels listed in eqs.(^); the 
difference is noticeable only for m cr{ ~ nih- The dotted curves show our estimates of the 
minimal signal that is visible on top of the smooth 77 background, for three different 
values of the integrated luminosity. The signal has been computed for m t = 150 GeV, 
tan/3 = 2, M 2 = 1.5m fl , m~ tL = 1.5m^ = 300 GeV, m P = 500 GeV and y, = -133 
GeV. 

Fig. 7 Cross section for production at the LHC. The solid line shows the total cross section 
multiplied with 0.01, and the dashed curves the 77 signal cross section after cuts, for 
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the two scenarios of figs. 2 and 4. The dotted curve shows the minimal cross section 
giving a significant signal after one year of nominal LHC operations, as defined in the 
text. 

Fig. 8 The dependence of the 77 signal cross section after cuts on the choice of scale Q in the 
gluon distribution functions, and on the QCD parameter A. The parameters are as in 
fig.5, with A = 1. 

Fig.9 The ratio J2n l-Rn(0)| 2 /|-Ri(0)| 2 , where n runs over all true stoponium bound states, 
defined by eq. (|TTD . This is a measure of the possible enhancement of the signal for 
stoponium production due to the production of excited states, as discussed in the text. 

Fig. 10 The region in the plane spanned by m a . and A t + fi cot/3 after one (solid) and five 
(dashed) years of running the LHC at full luminosity (£ = 100 fb _1 per year). The 
region in the top left corner is exlcuded by LEP searches for neutral Higgs bosons. 
The curves have been obtained for m t = 150 GeV, tan/? = 3, M2 = 1.67/77.^, m^ L = 
Tni R , /j, = 750 GeV and rap = 2 TeV, but depend little on these choices unless m a , ~ 
mp, as discussed in the text. 
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